1 and 0"^"^ heavy four-quark and molecule states in QCD 

R.M. Albuquerque^'' * F. Fanomezana*^ ** S. Narison'''', A. Rabemananjara'^ ** 

"Instituto de Fisica, Universidade de Sao Paulo, C P. 66318, 05389-970 Sao Paulo, SP, Brazil 
^ Laboratoire Particules et Univers de Montpellier, CNRS-IN2P3, Case 070, Place Eugene Bataillon, 34095 - MontpelUer, France. 
'^Institute of High-Energy Physics of Madagascar (iHEP-MAD), University of Antananarivo, Madagascar 



Abstract 

We estimate the masses of the 1 heavy four-quark and molecule states by combining exponential Laplace (LSR) and finite 
"energy (FESR) sum rules known perturbatively to lowest order (LO) in but including non-perturbative terms up to the complete 
dimension-six condensate contributions. This approach allows to fix more precisely the value of the QCD continuum threshold 
■(often taken ad hoc) at which the optimal result is extracted. We use double ratio of sum rules (DRSR) for determining the S U(3) 
breakings terms. We also study the effects of the heavy quark mass definitions on these LO results. The S U{3) mass-splittings of 
about (50 - 110) MeV and the ones of about (250 - 300) MeV between the lowest ground states and their 1st radial excitations 
are (almost) heavy-flavour independent. The mass predictions summarized in Table 4 are compared with the ones in the literature 
(when available) and with the three 7^(4260, 4360, 4660) and 7^(10890) 1 experimental candidates. We conclude (to this order 
approximation) that the lowest observed state cannot be a pure 1 four-quark nor a pure molecule but may result from their 
■mixings. We extend the above analyzes to the 0^^ four-quark and molecule states which are about (0.5-1) GeV heavier than the 
'"ZT' corresponding 1 states, while the splittings between the 0^^ lowest ground state and the 1st radial excitation is about (300-500) 
Q^MeV. We complete the analysis by estimating the decay constants of the 1 and 0^^ four-quark states which are tiny and which 
I exhibit a 1 /Mq behaviour. Our predictions can be further tested using some alternative non-perturbative approaches or/and at LHC^ 
some other hadron factories. 

D • 
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1. Introduction and a short review on the 1^"^ channel 

■ A large amount of exotic hadrons which differ from the "stan- 
cJard" cc chamonium and bb bottomium radial excitation states 
have been recently discovered in B-factories through J/if/n^n'' 
and Yn^n^ processes and have stimulated different theoretical 
interpretations. Most of them have been assigned as four-quarks 
and/or molecule states [3]. In previous papers [1, 2], some 
of us have studied, using exponential QCD spectral sum rules 
(QSSR) [4] 2 and the double ratio of sum rules (DRSR) [7] ^ the 
nature of the X(3872) 1^^ states found by Belle [11] and con- 
_firmedby Babar [12], CDF [13] and DO [14]. If it is a (c^)(^ 
four-quark or D - D* molecule state, one finds for nic - 1 .23 
'GeV [1] 4; 

Xc = (3925 ± 127) MeV , (1) 

corresponding to a fc-value common solution of the exponential 
Laplace (LSR) and Finite Energy (FESR) sum rules; 

Vf7 = (4.15 + 0.03) GeV, (2) 
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while in the Zj-meson channel, using mt, 
finds [1]: 

Xb = (10144+104)MeV with Vf7 ■■ 



4.26 GeV, one 



(10.4+0.02) GeV,(3) 



where a similar result has been found in [15] using another 
choice of interpolating current. However, in the case of the 
Xc(3872), the previous two configurations are not favoured by 
its narrow hadronic width (< 2.3 MeV), which has lead some of 
us to propose that it could be, instead, a. A - J/tfr-type molecule 
[2] described by the current: 

4 = (j)\cA"rcma75q), (4) 

where /!„ is the colour matrix , while g and A are coupling and 
scale associated to an effective Van Der Vaals force. In this case, 
the narrow width of the X^ is mainly due to the extra-gluon ex- 
change which gives a suppression of the order or^ compared to 
the two former configurations, if one evaluates this width using 
vertex sum rules. The corresponding mass is slightly lower than 
the one in Eq. (1) [2]: 

" = 0.96 + 0.03 X;' = (3768 ± 127) MeV .(5) 



X 



mol 



which (within the errors) also agree with the data. By assuming 
that the mass of the radial excitation X'q x ^Jt^, one can also 
deduce the mass-splitting: 



K-Xc^ 225 MeV ^X'-Xb^ 256 MeV , 



(6) 
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which is much lower than the ones of ordinary charmonium and 
bottomium states: 

il/{lS) - ifrilS) ^ 590 « T(25) - T(15) ^ 560 MeV, (7) 

and suggests a completely different dynamics for these ex- 
otic states. Comparing the previous results with the observed 
Zi,(10610) and Z;,(10650) states whose quantum numbers have 
been assigned to be l"^"^, one can conclude that these observed 
states are heavier than the 1st radial excitation of the Xfc(10.14) 
expected from QSSR to lowest order in ckj [1]. 

2. QCD Analysis of the 1 and 0"*"* channels 

In the following, we extend the previous analysis to the case 
of the 1 and 0^^ channels and improve some existing analysis 
from QCD (spectral) sum rules in the 1"~ channel [18, 19]. The 
results will be compared with the experimental 1 candidate 
states: 

F(4260), 7(4360), F(4660) , 7^(10890) (8) 

seen by Babar [16] and Belle [17, 23] and which decay into 

Jjif/n'^n^ and fn^n^ around the T(55) mass. These states 
cannot be identified with standard cc charmonium and bb bot- 
tomium radial excitations and have been assigned in the litera- 
ture to be four-quark or molecule states or some threshold ef- 
fects. 

• QCD input parameters 

The QCD parameters which shall appear in the following 
analysis will be the charm and bottom quark masses rricb, 
the light quark masses m^ j, the light quark condensates 
{qq) and {ss), the gluon condensates (g^G^) = ig^G'^yGg^) 
and (g^G^) = {g^fabcG%GlpGl^), the mixed condensate 
(qgcrGq) = {qgo-''^(Aa/2)G'^^yq) and the four-quark condensate 
piqq}^, where p indicates the violation of the four-quark vac- 
uum saturation. Their values are given in Table 1 and we shall 
work with the rurming Ught quark parameters: 

OT (r) = 

(- log v?a) 

{qq)(T) = -/i^(-log VtA) 

{qgcrGq){T) = -Mlfi][-\og ^Rk)~"^^\ (9) 

where px = -(1/2)(11 - 2n/3) is the first coefficient of the jS 
function for n flavours; and fiq are renormalization group 
invariant light quark mass and condensate [25, 26]. 

• Interpolating currents 

We assume that the Y state is described either by the lowest 
dimension (without derivative terms) four-quark and molecule 
DsD* vector currents given in Tables 2 and 3. Urilike the 
case of baryons where both positive and parity states can cou- 
ple to the same operator [24], the situation is simpler here as 
the vector and axial- vector currents have a well-defined quan- 
tum numbers to which are associated the 1"" (resp. 1"^^) states 



Table 1: QCD input par amete rs. For the heavy quark masses, we use the range 
spaimed by the rurming MS mass mg(Mg) and the on-shell mass from QCD 
(spectral) sum rules compiled in pages 602 and 603 of the book in [5] and 
recently obtained in Ref. [28]. The values of A and flq have been obtained from 

a^iMr) = 0.325(8) [29] and from the running masses: (m„ + mrf)(2) = 7.9(3) 
MeV [31]. The original errors have been multiplied by 2 for a conservative 
estimate of the errors. 



Parameters 


Values 


Ref. 


A(nf = 4) 


(324 ± 15) MeV 


[29, 30, 32] 


A(n/ = 5) 


(194 ± 10) MeV 


[29, 30, 32] 


m. 


(0.114 ± 0.021) GeV 


[5,31,32] 


rric 


(1.26 ~ 1.47) GeV 


[5, 28,31-34] 


rrib 


(4.17 ~ 4.70) GeV 


[5, 28,31-33] 




(263 ± 7) MeV 


[5,31] 


K = {ss}/{uu} 


(0.74 ± 0.06) 


[9] 


Ml 


(0.8 ± 0.2) GeV^ 


[35? , 36] 




(7 ± 2) X 10-2 GeV^ 


[28, 29, 37^3] 




(8.3 ± 1.0) GeV^ X (or.G^) 


[28] 


p = {q<iqq)l{qqf 


(2±1) 


[29, 35, 37] 



for the transverse part and the 0"^^ (resp. 0"") states for the lon- 
gitudinal part. In the case of four-quark currents, we can have 
two-types of lowest derivative vector operators which can mix 
through the mixing parameter b ^ . Another possible mixing can 
occur through the renormalization of operators [26, 27] though 
this type of mixing will only induce an overall effect due to the 
anomalous dimension which wiU be relevant at higher order in 
OTj but will disappear in the ratio of sum rules used in this paper. 
For the molecule current, we choose the product of local biUn- 
ear current which has the quantum number of the corresponding 
meson state. In this sense, we have only an unique interpolat- 
ing current. Observed states can be a mixing of diflFerent states 
associated to each choice of operators and their selection can 
only be done through the analysis of their decays [2] but this is 
beyond the scope of this paper. 

• The two-point function in QCD 

The two-point functions of the Yq {Q = c, b) (assumed to be a 

1 ^ vector meson) is defined as: 

U^^q) ^ ij d^xe'''mT[nx)fHO)m 

= -n(')(^2)( ,v _ ^ U^yq^f-^ , (10) 

r r 

where are the interpolating vector currents given Tables 

2 and 3. We assume that the Y state is described either by 
the lowest dimension (without derivative terms) four-quark and 
molecule D^D* currents given in Tables 2 and 3. The two in- 
variants, n*'' and n*^"', appearing in Eq. (10) are independent 
and have respectively the quantum numbers of the spin 1 and 
mesons. We can extract H^^ and n^'''(^^) or the corresponding 



""The 1++ four-quark state described by the axial-vector current has been 
analyzed in [1, 2]. 



spectral functions from the complete expression of n'^iq) by 
applying respectively to it the projectors: 



and 



q^qy 



(11) 



Due to its analyticity, the correlation function, n*'''^'(^^) in 
Eq. (10), obeys the dispersion relation: 



^ J Ami 



Im n^^'^s) 

as T — — + ■ 

s — q'^ — le 



(12) 



where Im n''''''(5) are the spectral functions. The QCD expres- 
sions of these spectral functions are given in Tables 2 and 3. 
1 Iq'^ terms discussed in [44, 45],which are dual to higher order 
terms of the QCD series will not be included here as we work 
to leading order. 



3. 1 four-quark state mass Fg, from QSSR 

In the following, we shall estimate the mass of the 1 four- 
quark state (Qq)(Qq) (Q = c, b and q = u, d quarks), hereafter 
denoted by Fq^. In so doing, we shall use the ratios of the 
Laplace (exponential) sum rule: 

??^^/(t) = -^Log J'' dte-''Umn^'\t) , (13) 



and ofFESR: 



rnFESR 



j^'' t/f f"ilmn(i>(f) 
£/f f«-iilmn(i)(f) 



n = 1 , 



(14) 



where f< is the hadronic (quark) threshold. Within the usual 
duality ansatz "one resonance" + 9{t - tc) x QCD continuum 
parametrization of the spectral function, the previous ratios of 
sum rules give: 



tnFESR 



(15) 



a) 



b) 
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m,=1.26 GeV 

-Vf7= 5.0+ 2mj 
■-/r=5.1+2mj 
■■Vi|=5.2+2mj 
*-Va~= 6.0+2mrf 

- \n7= 7.0+2mj 
<^Vir=7.5+2ni,, 







m,=1.47 GeV 

— -/rr=5.5+2m,; 
-*-^ft^=5.(i+2m,i 
■-- V^=5.8+2«i,, 

-- v7~= 6.5+2mj 
^Vt~=7-5+2n;,; 
- ^[^^^ 8.0+2m,; 
-^y/T^ 8.5+2m,, 





c) 




LSR 1.26 
FESR 1.26 
LSR 1.47 
FESR 1 .47 



For a discussion more closed to the existing literature which we 
shall test the reliability in the following, we start to work with 
the current corresponding to = 0. We shall discuss the more 
general choice of current when b is a free parameter at the end 
of this section. 

• The Ycd mass from LSR and FESR for the case b—0 

Using the QCD inputs in Table 1 , we show the r-behaviour of 
My^.j from 'R^'^ in Fig. la for - 1.26 GeV and for different 
values of tc- One can notice from Fig. la that the r-stability is 
obtained from V^c > 5.1 GeV, while the ff-stability is reached 
for Vf^ - 7 GeV. The most conservative prediction from the 
LSR is obtained in this range of fc-values for nic = 1 .26 GeV 
and gives in units of GeV: 

4.79 < Mf„, < 5.73 for 5.02 < < 7 and m,. = 1.26, 
5.29 < My„, < 6.11 for 5.5 < V^J < 7 andm,. = 1.47. (16) 



5.0 5.5 6.0 6.5 7.0 7.5 8.0 

VfT [GeV] 



Figure 1 : a) r-behaviour of My^^ (1 ) from *?^^^ for the cuneiit mixing parameter b = 0, 
for different values of and for - 1 .26 GeV; b) The same as a) but for = 1 .47 GeV; 
c) -behaviour of the LSR results obtained at the r-stability points and comparison with 
the ones from K™^* form,. = 1.26 and 1.47 GeV. 



We compare in Fig. lb), the ff-behaviour of the LSR results 
obtained at the r-stability points with the ones from R^^^'^ for 
the charm quark mass mc=1.23 GeV (running) and 1.47 GeV 
(on-shell). One can deduce the common solution in units of 
GeV: 

My^, = 4.814 for -y/fJ = 5.04(5) and m,. = 1.26, 

= 5.409 for Vf^ = 5.6 and m,. = 1.47 . (17) 

In order to fix the values of My^^ obtained at this lowest order 
PT calculations, we can also refer to the predictions of the J/ if/ 
mass using the LSR at the same lowest order PT calculations 
and including the condensate contributions up to dimension- 
six. We observe that the on-shell c-quark mass value tends to 
overestimate Mjj^ [2, 28]. The same feature happens for the 
evaluation of the X{\^^) four-quark state mass [1]. Though this 
observation may not be rigorous as the strength of the radia- 
tive corrections is channel dependent, we are tempted to take 
as a final result in this paper the prediction obtained by using 
the running mass mdmc) - 1262(17) MeV within which it is 
known, from different examples in the literature, that the PT se- 
ries converge faster [28] ^. Including different sources of errors, 
we deduce in MeV ': 

My,, = 4814(50)ai4),„,(2)A(17)r„,(2)c2(4)M2(13)G3(6)p 



^We plan to check this conjecture in a future publication when PT radiative 
corrections are included. 

'We consider this result as an improvement (smaller error) of the one e.g. 
in [18, 19] where only exponential sum rales have been used. However, the 
present error and the existing ones in the Hterature may have been underesti- 
mated due to the non-inclusion of the unknown PT radiative coiTections and 
some eventual systematics of the approach. 
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= 4814(57). 



(18) 



Using the fact that the 1st FESR moment gives a correlation 
between the mass of the lowest ground state and the onset of 
continuum threshold tc, where its value coincide approximately 
with the value of the 1st radial excitation mass (see e.g. ref. 
[39] and some other examples in [5]), we shall approximately 
identify its value with the one of the radial excitation. In or- 
der to take into account the systematics of the approach and 
some eventual small local duality violation advocated by [46] 
which can only be detectable in a high-precision analysis like 
the extraction of as from r-decay [29, 47], we have allowed fc 
to move around this intersection point. Assuming that the mass 
of the radial excitation is approximately ^/t^, one can deduce 
the mass-splitting: 



M'y -My 



226 MeV 



(19) 



which is similar to the one obtained for the X(l^^) four-quark 
state [ 1 ] . This splitting is much lower than the one intuitively 
used in the current hterature: 

M^(25)-M^(15)^ 590MeV, (20) 

for fixing the arbitrary value of tc entering in different Borel 
(exponential) sum rules of the four-quark and molecule states. 
This difference may signal some new dynamics for the exotic 
states compared with the usual cc charmonium states and need 
to be tested from some other approaches such as potential mod- 
els, heavy quark symmetry, AdS/QCD and lattice calculations. 



a) 
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T [GeV--] 



nii=4.17 GeV 

— Vi"- 11.0+2m^ 
*Vi~= 11.2+2mj 
■■■^ft7= 11.5+2hi^ 
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I, = 14.5+2m,, 
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Figure 2: a) r-behaviour of My^^(1 ) from 'R^j'^ for the cuirent mixing parameter h = 0, 
for different values of t^- and for m;, - 4. 17 GeV. b) /^.-beiiaviour of tlie LSR results obtained 
at the T-stability points and comparison with the ones from fi^jf^^ for m/, = 4.17 and 4.70 
GeV. 



• The Ybd mass from LSR and FESR for the case b — 
Using similar analysis for the /7-quark, we show the r-behaviour 
of '??^^*(t) inFig. 2a for = 4.17GeV and for different values 
of tc- In Fig. 2b, the same analysis is shown for m/, = 4.70 GeV. 
The most conservative result from the LSR is (in units of GeV) 
is: 

11.0 < My,, < 12.4 for 11.2 < < 14.5 and = 4.17, 

12.1 < My, , < 13.4 for 12.2 < -y^ < 15.5 and mt = 4.70, 



where the lower (resp. higher) values of tc correspond to the 
beginning of r (resp. fc)-stability. We compare in Fig. 2b), 
the fc-behaviour of the LSR results obtained at the r-stability 
points with the ones from ^^J'^* for the b quark mass mt-A-Al 
GeV (running) and 4.70 GeV (on-shell). One can deduce the 
common solution in units of GeV: 



My,, = 11.26 for ^f^c■ 
= 12.09 for ^f^c■■ 



11.57(7) and mi =4.17 
12.2 and OTA = 4.70. (21) 



One can notice, like in the case of the charm quark that the value 
of the on-shell quark mass tends to give a higher value of My,, 
within this lowest order PT calculations. Considering, like in 
the case of charm, as a final estimate the one from the running b- 
quark mass nihimb) - 4177(1 1) MeV [28], we deduce in MeV: 

My,, = 11247(45),,.(8)„,(2)A(15)r,„(l)G^(l)M2(l)G3(5)p 
= 11256(49). (22) 

From the previous result, one can deduce the approximate value 
of the mass-splitting between the 1st radial excitation and the 
lowest mass ground state: 



M'y -M- 



250 MeV 



(23) 



which are (almost) heavy-flavour independent and also smaller 
than the one of the bottomium splitting: 



My(25) - Mt(15) ^ 560 MeV 



Effect of the current mixing b on the mass 



(24) 



In the following, we shall let the current mixing parameter 
b defined in Table 2 free and study its effect on the results ob- 
tained in Eqs. (18) and (22). In so doing, we fix the values of t 
around the r-stability point and tc around the intersection point 
of the LSR and FESR. The results of the analysis are shown in 
Fig. 3. We notice that the results are optimal at the value b - Q 
which a posteriori justifies the results obtained previously for 
Z7 =0. 



a) 



b) 
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Figure 3 : a) Zj-behaviour of My^,^ for given values of r and tc and for i 
the same as a) but for ^^y,, and for fli^ =4.17 GeV. 



1.26 GeV; b) 
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• Effect of the current mixing b on the decay constant fygj 

For completing the analysis of the effect of b, we also study 
the decay constant fy^d defined as: 



(25) 



We show the analysis in Fig. 4 giving My^j and the correspond- 
ing fc obtained above. One can deduce the optimal values at 
b^O: 



fy^, ^ 0.08 MeV and /y,, ^ 0.03 MeV 



(26) 



which are much smaller than f„ = 132 MeV, ^ ^ 215 MeV and 
fo - /b=203 MeV [48]. On can also note that the decay con- 
stant decreases like 1 /Mq which can be tested in HQET or/and 
lattice QCD. 
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Figure 4: a) r-behaviour of fy^,^ for given values oi h = and l^- and for = 1.26 GeV; 
b) the same as a) but for fy^^ and for m/, - 4.17 GeV; c) /j-beiiaviour of fy^^ for given 
values of r at the stability and d) the same as c) but for fy. . 



• S U(3) breaking for My^^ from DRSR 

We study the ratio My^JMy^^ using double ratio of LSR 
(DRSR): 



sd 



4 



LSR 



where Q = c, b . 



(27) 



We show the r-behaviour of rj^ and rj^ respectively in Figs. 5a 
and 5b for m^ = 1 .26 GeV and nih -4.17 GeV for different val- 
ues of tc- We show, in Fig. 5c and Fig. 5d, the fc-behaviour of 
the stabilities or inflexion points for two diff'erent values (run- 
ning and on-shell) of the quark masses. One can see in these 
figures that the DRSR is very stable versus the t^ variations in 
the case of the running heavy quark masses. We deduce the 
corresponding DRSR: 



= 1.018(l)„,(5)„,(2).(2)s„(l)p, 

= 1.007(0.5)„,(2)„X0.5).(l)s„(0.3)p, 



(28) 



respectively for = 5.1 and 11.6 GeV. Using the results for 
1.06 R 

a) 
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Figure 5: a) r-behaviour of r^^ for the current mixing parameter b = Q, for different 
values of tc and for nic = 1.26 GeV; b) r-behaviour of r'^^ for different values of tc and for 
mi, —4.17 GeV; c) /^.-behaviour of the inflexion points {or minimas) of r^^^ from Fig 5a; d) 
The same for the b quark using r^^ from Fig 5b. 



Yqd in Eqs. (18) and (22) and the values of the S U(3) breaking 
ratio in Eq. (28), we can deduce the mass of the Yq, state in 
MeV: 



My„ = 4900(67) , My,, = 1 1334(55) , 
leading to the S U{3) mass-splitting: 

AMf^ ^ 87 MeV ^ AM^^ ^ 78 MeV , 
which is also (almost) heavy-flavour independent. 



(29) 



(30) 



4. 1 molecule masses from QSSR 

• The D*j^^^Dd{s) and B*j^^^Bj(s) molecules^ 

Like in the previous case, we use LSR and FESR for studying 
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Figure 6: a) r-behaviour of M^* for diiferent values of and for mc = 1.26 GeV; b) 
T-behaviour of M^'^^ for different values of and for nij, =4.17 GeV; c) f^-behaviour of 
the extremas in r of and for = 1.26- 1.47 GeV; d) the same as c) but for Mg^e^ 

and for mt = 4.17 - 4.70 GeV. 



the masses of the D*J)d and B'^Bd and DRSR for studying the 
S U (3) breaking ratios: 



Mm, 



M, 



B ^ ^B-B, 
''sd - 



(31) 



BjBd 



We show their r-behaviour for different values of tc and for 
nic - 1.26 GeV and mt, = 4.17 GeV respectively in Figs. 6a,b 
and 7a,b. The f^-behaviour of the r-minimas is shown in Fig. 
6c,d for the masses and in Fig. 7c,d for the S U{3) breaking 
ratios. Using the sets (nic - 1.26 GeV, -\/f7 = 5.58 GeV) and 
(mb - 4.17 GeV, ^/t^ - 11.64(3) GeV) common solutions of 
LSR and FESR, one can deduce in MeV: 

Mo-D, = 5268(14)„J3)A(19)r,„(0)c=(0)M2(2)c3(5)p, 
= 5268(24) , 

Ms'^s, = 1 1302(20)^9)™, (2)a(19)s„(0)g2(0)m^(1)cK5)p 



Hereafter, for simplifying notations, D and B denote tlie scalar and 
mesons. 
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Figure 7: a) r-behaviour of for different values of f,. and for m^- - 1.26 GeV; b) r- 
behaviour of f^^ for different values of and for m/, - 4.17 GeV; c) /^.-behaviour of the 
inflexion points (or minimas) of from Fig 7a; d) the same for the b quark using r^^ from 
Fig 7b. 



11302(30), 

1.018(l),„/4)m,(0.8),(0.5)s„(0.2)p(0.1)G3 , 
rf, = 1.006(l)„,„(2)m,(l).(0.5)s„(0.2)p(0.1)G3 . (32) 
Using the previous results in Eq. (32), one obtains in MeV : 

Md;d, = 5363(33) , Mb-^b, = 1 1370(40) , (33) 
corresponding to a 5 f/(3) mass-splitting: 



• sd 



AM 



sd 



95 MeV X AM' 



sd 



68 MeV 



(34) 



These results for Mdd' are in the upper part of the range given 
in [18] due both to the smaller values of - 1.23 GeV and 
^/t^-5.5 GeV used in that paper Though the DD* molecule 
mass is above the DD* threshold which is similar to the e.g. 
the case of the nn continuum and p-meson resonance in e^e' 
to the 1=1 hadrons channel, one expects that at the r-stability 
point or inside the sum rule window, where the QCD continuum 
contribution is minimum while the OPE is still convergent, the 
lowest ground state dominates the sum rule. 

• The J lips 2 and TS2 molecules 

Combining LSR and FESR, we consider the mass of the J/1//S2 
and TS2 molecules in a colour singlet combination, where ^2 = 



a) 
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c) 
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Figure 8: a) r-behaviour of Mj/^s^ for different values of and for = 1.26 GeV; b) 
T-beiiaviour of Mrs 2 for ditferent values of and for m/, - 4. 17 GeV; c) -behaviour of 
the extremas in r of Mj/^s^ for = 1.26 - 1.47 GeV; d) the same as c) but for Mrs^ for 
=4.17 -4.70 GeV. 



Figure 9: a) r-behaviour of for different values of and for = 1.26 GeV; b) 
T-behaviour of r^^ for different values of and for - 4.17 GeV; c) rc-behaviour of 
the extremas in t of r^^ for m^- = 1.26 - 1.47 GeV; d) the same as c) but for r^^ for 
Hij =4.17 -4.70 GeV."' 



MM + dd is a scalar meson ^. In so doing, we work with the LO 
QCD expression obtained in [19]. We show the results versus 
the LSR variable t in Fig. 8a,b. The fc-behaviour of different 
T-extremas is given in Figs. 8c,d from which we can deduce for 
the running quark masses for ^/t^ = 5.30(2) and 10.23(3) GeV 
in units of MeV: 

Mj/^s, = 5002(20)^.(8),„.(2)A(19)f,„(9)c2(0)M2(0)c3(6)p 

= 5002(31), 

Mrs, = 10015(20),.(9)„,(2)A(16)a„(17)c2(0)^2(0)c3(5)p 

= 10015(33) . (35) 

The splitting (in units of MeV) with the first radial excitation 
approximately given by ^/t^ is: 



'sd 



Mrs, 
Mrs 2 



= 1.011(1)„,(2)„,(0.2),, 



(37) 



where 5 3 = is a scalar meson. Then, we obtain in MeV: 

Mj/^s, = 5112(41) , Mrs, = 10125(40) , (38) 
corresponding to the S U(3) mass-splittings: 



AM' 



y/iA 

sd 



AM 



sd 



110 MeV , 



(39) 



M' 



JlfS2 



M 



J/^S2 



298 



M' 



rsi 



Mrs, ~ 213. (36) 



In the same way, we show in Figs. 9 the r and tc behaviours of 
the S U{3) breaking ratios, from which, we can deduce: 



Mj/^Si 



'sd 



M 



J/lf,S2 



= 1.022(0.2),„^.(5)„,,(2),, 



'The low-mass ;r+;r" invariant mass due to the cr meson is expected to result 
mainly from its gluon rather than from its quark component [49, 50] such that 
an eventual quark-gluon hybrid meson nature of the is also possible. 



The mass-splittings in Eq. (39) are comparable with the ones 
obtained previously. 

Doing the same exercise for the octet current, we deduce the 
results in Table 4 where the molecule associated to the octet 
current is 100 (resp. 250) MeV above the one of the singlet 
current for J/i// (resp T) contrary to the 1^^ case discussed in 
[2]. The ratio of S U(3) breakings are respectively 1.022(5) and 
1.010(2) in the c and b channels which are comparable with 
the ones in Eq. 37. When comparing our results with the ones 
in Ref. [19], we notice that the low central value of Mjj^s, 
obtained there (which we reproduce) corresponds to a smaller 
value of nic - 1 .23 GeV and mainly to a low value of ^Jt^ = 5.1 
GeV which does not coincide with the common solution ^[^c = 
5.3 GeV from LSR and FESR. On the opposite, the large value 
of Mrs, - 10.74 (resp. 1 1 .09) GeV obtained there corresponds 
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to a too high value ^|t^ =11.3 (resp. 11.7) GeV compared with 
the LSR and FESR solution = 10.23 (resp. 10.48) GeV for 
the singlet (resp. octet) current. 



a) 
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- 



b) 



c) 



d) 
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Figure 10: a) r-behaviour of for the cuirent mixing parameter /? - 0, for different 
value.s of and for m^. = 1.26 GeV ; b) r-behaviour of 1^^ for different values of and 
for m^j = 4.11 GeV; c) ^^.-behaviour of the extremas in r of for = 1.26 - 1.47 GeV; 
d) the same as c) but for Y?. for MJt = 4. 17 - 4.70 GeV. 
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Figure 11: a) /j-behaviour of MyO for given values of r and r^. and for nic 
b) the same as a) but for and for m/, = 4. 17 GeV. 



1.26 GeV; 



One can notice that the splittings between the lowest ground 
state and the 1st radial excitation approximately given by 
is in MeV: 



M'^ - Myo^ ^ 375 , M'^ - Myo^ ^ 464 



(41) 



which is larger than the ones of the 1 states, comparable with 
the ones of the J/i// and T, and are (almost) heavy-flavour in- 
dependent. We show in Fig 11 the effect of the choice of b 
operator mixing parameter on the mass predictions, indicating 
an optimal value at b - 0. For completeness, we show in Fig. 
12 the T and b behaviours of the decay constants from which 
we deduce: 



/yo 0.12 MeV and /yo ^ 0.03 MeV , 



(42) 



which are comparable with the ones of the spin 1 case in Eq. 
(26). 



5. O*"*" four-quark and molecule masses from QSSR 

In the following, we extend the previous analysis to the case of 
the 0"^"^ mesons. 

• mass and decay constant from LSR and FESR 

We do the analysis of the F"^ and Y^^ masses using LSR and 
FESR. We show the results in Figs 10 for the current mixing 
parameter b - Q from which we deduce in MeV, for the running 
quark masses, and respectively for - 6.5 and 13.0 GeV 
where LSR and FESR match: 

Myo = 6125(16),„,.(7)A(44)fl„(12)G:(14)p 

cd 

= 6125(51) MeV, 

Myo^ = 12542(22),X13);„,(l)A(7)r,„(34)G2(2)p 

= 12542(43) MeV. (40) 



• S U(3) breaking for from DRSR 

We show in Figs. 13 the t and tc behaviours of the S U{3) break- 
ing ratios for the current mixing parameter b -Q: 

/f.^: Q.cb, (43) 

Qd 

from which we deduce: 

= 1.011(2)„X3.8)m,(1.4),(l)r»(0.7)p, 
rf, = 1.004(l)„X1.7)m,(0.3),, (44) 

leading (in units of MeV) to: 

Myo = 6192(59), Myo = 12592(50) , (45) 

and the S U(3) mass-splittings: 

AM^. ^67 ~ AM^' =^ 50 MeV . (46) 
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Figure 12: a) r-behaviour of f^o for given values of /j = and and for = 1.26 

cd 

GeV; b) the same as a) but for fyo and for - 4.17 GeV; c) /7-behaviour of fyo at the 
T-stability and for a given value of d) the same as c) but for /j,o 



Figure 13: a) r-behaviour of r^^ for the current mixing parameter = 0, for different 
values of t^. and for rUc = 1 .26 GeV; b) r-behaviour of r^^ for different values of fc and for 
nib - 4.17 GeV; c) ft-behaviour of the extremas in r of r^^ for - 1.26 - 1.47 GeV; d) 
the same as c) but for f™' for = 4.17 - 4.70 GeV. 



• MojDj ond MbjB,, from LSR and FESR 

We show the t and behaviours of the masses M^jOd ^nd 
MbjBj in Figs 14. Like in previous sections, we consider as 
a final result (in units of MeV) the one corresponding to the 
running masses for ^/t^. - 6.25(3) and 12.02 GeV: 

Mdm = 5955(24)^.(14)^^(5)a(36)s„(4)c2(4)c3(12)p 

= 5955(48) , 

Mb,b, = 11750(12),„,X4)A(35)f,„(7)G2(3)G3(12)p 

= 11750(40) (47) 

One can notice that the splittings between the lowest ground 
state and the 1st radial excitation approximately given by 
is in MeV: 

M'o^Od - Mdm 290 , M',^^,^^-MB,Bd^21Q, (48) 

which, like in the case of the 1 states are smaller than the ones 
of the J 1 4/ and Y, and almost heavy-flavour independent. 

• S U(3) breaking for M^ ^^ and M^ b, from DRSR 

We show in Fig. 15 the r- behaviour of the SU(3) mass ra- 
tios for different values of tc and the tc behaviour of their t- 



extremas. Therefore, we deduce: 

C = = 1.015(l)„/4)m,(2),(l)s„(0.5)p, 

^DdDd 

rfd = = 1.008(l),„^.(4)m,(2),(l)s„(0.5)p. (49) 

iMB.Bd 

Using the previous values of MojOj and MbjBj, we deduce in 
MeV: 

Md,d, = 6044(56) , Mb,b, = 1 1844(50) , (50) 
which corresponds to a SU{3) splitting: 

AM^f ^ 89 MeV ^ AMff ^ 94 MeV . (51) 

6. Summary and conclusions 

We have studied the spectra of the 1 and 0^^ four-quarks 
and molecules states by combining Laplace (LSR) and finite en- 
ergy (FESR) sum rules. The S U(3) mass-splittings have been 
obtained using double ratios of sum rules (DRSR). We con- 
sider the present results as improvement of the existing ones in 
the literature extracted only from LSR where the criterion for 
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Figure 14: a) r-behaviour of Mo^Dj for different values of /j and for = 1 .26 GeV; b) 
T-behaviour of Mb^^b^ for different values of tc and for nih - 4. 17 GeV; c) rc-behaviour of 
the extremas in r of ^d^d^ 'ind for = 1 .26 — 1 .47 GeV; d) the same as c) but for Mb^b^ 
and for = 4. 17 - 4.70 GeV. 



Figure 15: a) r-behaviour of rf^ for different values of and for - 1.26 GeV; b) 
T-behaviour of r^^j for different values of and for mi, = 4.17 GeV; c) ?c-behaviour of the 
inflexion points {or minimas) of from Fig 7a; d)the same for the b quark using r^^ from 
Fig 7b. 



fixing the value of the continuum thresholds are often ad hoc 
or based on the ones of the standard charmonium/bottomium 
systems mass-splittings which are not confirmed by the present 
analysis. Our results are summarized in Table 4. We find that : 

• The three 7^(4260, 4360, 4660) 1"" experimental candi- 
dates are too low for being pure four-quark or/and molecule 
DD* and Jlij/So states but can result from their mixings. The 
7^(10890) is lower than the predicted values of the four-quark 
and BB* molecule masses but heavier than the predicted fS2 
and "XSt, molecule states. Our results may indicate that some 
other natures (hybrids, threshold effects,...) of these states are 
not excluded. On can notice that our predictions for the masses 
are above the corresponding meson-meson thresholds indicat- 
ing that these exotic states can be weakly bounded. 

• For the 1 , there is a regularity of about (250-300) MeV for 
the value of the mass-splittings between the lowest ground state 
and the 1 st radial excitation roughly approximated by the value 
of the continuum threshold yft^ at which the LSR and FESR 
match. These mass-splittings are (almost) flavour-independent 
and are much smaller than the ones of 500 MeV of ordinary 
charmonium and bottomium states and do not support some ad 
hoc choice used in the literature for fixing the -values when 
extracting the optimal results from the LSR. 

• There is also a regularity of about 50-90 MeV for the S U(3) 



mass-splittings of the different states which are also (almost) 
flavour-independent. 

• The spin states are much more heavier (> 400 MeV) than 
the spin 1 states, like in the case of hybrid states [5]. 

• The decay constants of the 1 and 0^^ four-quark states ob- 
tained in Eqs (26)and (42) are much smaller than f„, fp and 
fD,B- Unlike /g expected to behave as 1/ -\JMq, the four-quark 
states decay constants exhibit a 1 /Mq behaviour which can be 
tested using HQET or/and lattice QCD. 

It is likely that some other non-perturbative approaches such 
as potential models, HQET, AdS/QCD and lattice calculations 
check the previous new features and values on mass-splittings, 
mass and decay constants derived in this paper We also expect 
that present and future experiments (LHCb, Belle, Babar,...) 
can test our predictions. 
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Table 4: Masses of the four-quark and molecule states from the present analy- 
sis combining Laplace (LSR) and Finite Energy (pESR). We have used double 
ratios (DRSR) of sum rules for extracting the S U(3) mass-splittings. The re- 
sults correspond to the value of the running heavy quark masses but the S f7(3) 
mass-splittings are less affected by such definitions. As already mentioned in 
the text for simplifying notations, D and B denote the scalar and Bq mesons. 
The errors do not take into account the unknown ones from PT corrections. 



States 




States 




Four-quarks 


1- 




0*+ 


Ycd 


4818(27) 


cd 


6125(51) 


Y 

^ cs 


4900(67) 


yO 

CS 


6192(59) 


Yhd 


11256(49) 


^bd 


12542(43) 


Yi„ 


11334(55) 


as 


12592(50) 


Molecules 


1- 






D-Pd 


5268(24) 


DdD, 


5955(48) 


DID, 


5363(33) 


D,D, 


6044(56) 




11302(30) 


BdBd 


11750(40) 




11370(40) 


BsBs 


11844(50) 


Singlet current 


1- 


Octet current 


1- 




5002(31) 




5118(29) 




5112(41) 




5231(40) 




10015(33) 




10268(28) 




10125(40) 




10371(45) 
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Table 2: QCD expression of the Four-Quark Spectral Functions to lowest order in ccg and up to dimension-six condensates: Q = c, i> is the heavy quark field. 
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Table 3: QCD expression of the Molecule Spectral Functions to lowest order in Os and up to dimension-six condensates: Q = c,bis the heavy quark field, while g' 
and A' are coupling and scale associated to an effective Van Der Vaals force. 
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